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ABSTRACT 

We investigate the one-dimensional interaction of a relativistic jet and an external medium. Rela- 
tivistic magnetohydrodynamic simulations show an anomalous boost of the jet fluid in the boundary 
layer, as previously reported. We describe the boost mechanism using an ideal relativistic fluid and 
magnetohydrodynamic theory. The kinetic model is also examined for further understanding. Simple 
scaling laws for the maximum Lorentz factor are derived, and verified by the simulations. 
Subject headings: galaxies: jets — magnetohydrodynamics (MHD) — methods: numerical — rela- 
tivistic processe — shock wavess 



1. INTRODUCTION 

Relativistic jets are considered in various contexts in 
high-ener gy astrophysics, su ch as active galact ic nuclei 
(AGNs) jUrry fc Padovan i||1995| lFerrari| |1998| , micro- 
quasars ^Mirabel & Rodriguez 1999), and potentiall y 
gamma-ray bursts (GRBs) ( ,Piran,2005, |Meszaros||2006[ ). 
The interaction between fast movmg jets (the relevant 
Lorentz factors are 7jet ~ 10-20 in AGNs and ^jet ^ 10^ 
in GRBs) and the surrounding medium is very impor- 
tant to understand global dynamics of the jet system, 
because it is related to the mass, momentum, and energy 
transport across the boundary layers. In this context, 
develop ment of velocity shear in stabilities h as been of 
interes t d Turland fc S cheuer| ([1976 ) ; Blandford fc Pringle 
11976); Ferrari ct al. (1980); Birkinshaw (1991); Bodo et 



|al., (2004) ; ^Osmanov et al.^ (2008) and references therein] . 
Moreover, a relativistic jet-medium boundary is a po- 
tential site of high energy particle accelerat ion as well 
( [Ostrowski 2000; Stawarz fc Ostrowski||2002[ )■ 
Recently, it has been reported that the jet-medium 
interaction is more complex than thought even in the 
simplest one-dimensional (ID) case. Raising a Riemann 
problem of rela tivistic hydrodynamics (RHD), Aloy fc 



Rezzolla (2006) showed that the tangential hydro dy- 



namic velocity and the relevant Lorentz factor (7_bl) in 
the boundary layer are anomalousl y accelerated {ibl > 
^jet) when the jet is over-pressured. Mizuno et al. (20081 
studied relativistic magnetohydrodynamic (RMHD) ef- 
fect, and reported that th e perpendicular magnetic field 
enhances the boost effect. 
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Komissarov et al. 
cussed a similar tangential boost in their R. 
lation of the coUapser jet. Such anomalous boost effect 
may be responsible for incr easing the jet's Lor entz fac- 
tor \k\ay fc Rezzol la 2006; Mizuno et aI1|2008D and for 
modulating t he raHiativc signature of the jet (Aloy fc 
Mimica|2008 ) . However, its physical mechanism remains 
unclear, and therefore no quantitative analysis has been 
performed. 

In this paper, we study the mechanism of the anoma- 
lous boost by using RHD/RMHD simulations and an 
analytic theory. In Section 2, we describe the problem 
setup. In Section 3, we present the simulation results. In 
Section 4, we construct an RHD/RMHD theory of the 
problem. In Section 5, we additionally discuss kinetic 



aspects. The last section Section 6 contains discussions 
and summary. 

2. PROBLEM SETUP 



Following earlier works ( Aloy fc Rezzolla|2006 Mizuno 
et al.||2008| , we study a ID Riemann problem in a jet- 
like configuration, which is schematically illustrated in 
Figure [T] A jet travels upward in the -|-z-direction in a 
stationary ambient medium. An interaction between the 
jet and the medium is considered in the x-direction, and 
we assume dy = dz =0. Initially they are separated by 
a discontinuity and we study the time evolution of this 
ID system. 



We employ the following ideal RMHD equations ( Anile 
1989 1. For convenience we set c = 1 and employ Lorentz- 



Heaviside units such that all (47r)^/^ factors disappear. 

dt{ip) + V ■{^pv) = Q (1) 
dtm + V • {-i'^wtvv bb+ pj) = (2) 
dt£ + V ■m = Q (3) 
dtB + V ■ {vB - Bv) = (4) 
E + vxB = (5) 



{E X B) 
B)v 

TpJ{T -I) + 9 



(6) 



m — y^wtV — bob = phv 
£ ^j'^wt - bobo -~pt 
b ^{B/j)+-f{v 
b„^j{v-B) 

Wt — ph + b^ = p -\- 

Pt =Pg + 56' 

In the above equations, 7 is the Lorentz factor, p is the 
proper mass density, v is the velocity, m is the momen- 
tum density, £ is the energy density, Wt is the total en- 
thalpy, h is the specific enthalpy, pt is the total pressure, 
Pg is the gas pressure, T is the gas temperature including 
the Boltzmann constant, and b" = (60, b) is the covariant 
magnetic field. Note that 6^ = b°'ba = B^/j^ + {v-B)^ = 
{B^ — E"^) is a Lorentz invariant. We use an equation of 
state with a constant polytropic index of F = 4/3. 

We developed an RMHD code to numerically solve 
the problem. W e employ a relativistic HLLD scheme 
([Mignone et~al][2009_; Miyoshi fc Kusano] [2005]) , which 
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Figure 1. Our jet geometry. We consider the jet in the left side 
(L) and the ambient medium in the right side (ij). 



considers multiple states inside the Riemann fan in or- 
der to resolve discontinuities better. We interp olate the 
spati al profile by a monotonized central limiter ( van Leer 
1977 1 and solve the temporal evolution by the second 
order total variation diminishing (TVD) Runge-Kutta 
method. Relativistic primitive variables are recovered 



by 



Mignone & McKinney (2007rs inversion scheme, 
he model parameters are presented in Tabic [T] The 
subscripts L and R denote the properties in the two re- 
gions (L for the left side or the jet, and R for the right 
side or the ambient medium). The Lorentz factor of the 
jet is set to 7jet = 7. We initially set = 0. In our ID 
configuration this automatically means = all the 
time. This condition B^ = allows us to simplify the 
numerical scheme, because a five wave HLLD problem 



is redu ced to a three-wave problem (see [Mignone et al. 
(|2009[), Section 3.4.1). The first model HI has no mag- 



netic fields (RHD). The other two models contain mag- 
netic fields inside the jet: the jet-aligned magnetic field 
{B^: model Ml) and the out-of-plane magnetic field {By-. 
model M2). Importantly, the total pressure pt^L is set to 
the same. These RMHD models are analog ous to the 
"poloidal " (Ml) and "toroidal" (M2) cases in |Mizuno et 
at] ( |2008[ ). The spatial domain of -0.2 < a: < U.2 is re- 
solved by 6400 grids. A ll simulation results are checke d 
by an analytic solver by Giacomazzo & RezzoUa (2006). 



3. RESULTS 

Shown in Figure [2]3 are simulation results of the model 
HI at t = 0.2. One can recognize a three-wave struc- 
ture: (1) a leftward rarefaction wave {x ^ —0.02), (2) 
a contact discontinuity that separates the jet and the 
ambient medium {x ~ 0.02), and (3) a right-going for- 
ward shock (x ~ 0.12). The system exhibits a self-similar 
evolution as those waves propagate in time. Numerical 
errors are negligible, thanks to the high resolution and 
the stable numerical scheme. In the rarefaction region 
between (1) and (2), the Lorentz factor of the fluid grad- 
ually increases from 7jet = 7, and then it reaches to the 
maximum (^ 11.7) at the left vicinity of the contact dis- 
continuity. This is consistent with the anomalous boost 
demonstrated in previous works. Hereafter, we denote 
this boosted region as the "boundary layer" and define 
the relevant Lorentz factor in the flat region ^bl- The 
tangential velocity increases there, as shown in the small 
box in Figure [2]a. 

The RMHD models evolve similarly as the RHD model 
HI does. Figure [2b compares the pressure profiles of the 
three models, and Figure [2|c shows the profiles of the 
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Figure 2. (Color online) (a) Simulation result of model HI at 
t = 0.2. The fluid Lorentz factor 7, the gas pressure pg, the normal 
velocity Vx, and the tangential velocity are presented. The 
tangential velocity in the boosted region is also zoomed up in 
the small box. (b) The total pressure pt {solid lines) and the gas 
pressure pg {dashed lines) in models HI {black), Ml {red thin line), 
and M2 {blue thick line) at t = 0.2. (c) The Lorentz factor 7 in 
three models at i = 0.2. The small numbers indicate the Lorentz 
factors in the relevant flat regions. 



Lorentz factor. Since the jet contains the magnetic field 
in the RMHD cases, the rarefaction wave fronts propa- 
gate faster than the RHD case, because the Alfven speeds 
(^ c in the proper frames) are faster than the sound 
speed (cs ^ c/v^ in the proper frame). One can also see 
the tangential discontinuities between the jet and the am- 
bient medium {x ^ 0.03 in model Ml, x ~ 0.01 in M2), 
where the magnetic pressure disappears and the gas pres- 
sure suddenly increases to maintain the total pressure. 
The anomalous boost similarly takes place on the jet side 
of the those discontinuities. The forward shocks are just 
out of sight from fig ures in RMHD cases. As reported by 
[Mizuno et al. ( 2008 1 , the model M2 with a perpendicular 
magnetic fiela^^B^ exhibits stronger boost (7sl ^ 16.2) 
than the model Ml with a parallel magnetic field (i?z) 
(7SL - 11.1). 



4. ANALYTIC THEORY 

4.1. RHD theory 

In this section we study the mechanics of the anoma- 
lous boost problem. First we examine the RHD case. 
Combining the momentum equation (Equation [2]) and 
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Note. — Models and parameters for the Riemann problems. The subscript L denotes the jet (left side) properties and R for the ambient medium 
(right side). In addition, two parameter surveys are performed by changing = (10"'^, 10~^) and Pg,ij = (3,0.3,0.1). 



the e nergy equation (Equation p| (Sakai & Kawata 
T980|, 



dtij^phv) +v(y ■ {"f^phv)] + ^^ph{v ■ V)v + Vpg = 0, 



t>5t(7 ph) - vdtPg + t»( V • (7 phv) = 0, 



we obtain 



7 Ph 



Dv 

Dt 



dt 



(7) 



Since dz — 0, the anomalous boost obviously comes from 
the last term, ph{D/Dt)vz ~ —dtPg. This term has no 
Newtonian counterpart, it is certainly a relativistic effect. 
In usual contexts, the term slows down the fluid bulk 
acceleration in the high-temperature regime {pg > p), 
as if the relativistic pressure increases the inertia. In 
this case, since the pressure decreases in the rarefaction 
region, the force in the last term boosts the fluid in the 
z-direction, until the fluid element reaches the constant- 
pressure region. We see that the term coverts excess 
internal energy to the energy of the bulk motion. 

Next, we arrange the momentum equation (Equation 
[2]) in the following way. 

7p(9t + V ■ V){'yhv) + dt{^ip) + V • {-ipv) ^hv = -Vpg. 

Using Equation [1] we obtain 
D 



IPj^ilhvz) 



0. 



(8) 



Thus, the specific momentum (the momentum density 
per the gas density in this frame) remains constant as it 
should be. This is because no external forces accelerate 
the fluid, and because the ideal fluid assumption docs 
not allow momentum transport in its own frame. We 
confirmed that "fhvz is well conserved in both sides in 
the simulation. 

In model HI, the jet velocity is initially relativistic 
{vz,L ~ 1), and then we expect 



7/1 ^ const. 



(9) 



in the rarefaction region. The behavior of Equation |9] is 
controlled by the gas temperature, T = (pg/p). When 
the gas is cold (T <C 1), both the specific enthalpy h ^ 
1 and the Lorentz factor 7 remain constant; no boost 
occurs. When the gas is relativistically hot (T ^ I), 
h ^ 4T becomes a function of T. In this limit, we find 

7T = jipg/p) - const. (10) 

We see that the Lorentz factor increases when the rela- 
tivistic temperature decreases. Physically this is relevant 



to the temporal decrease of the pressure (Equation [?]). 
Combining with the polytropic law {pgP^^ = const.), we 
obtain the following relations. 



jp^ ^ ^ const. 
7^^'"""'^-'/'" const. 



(11) 
(12) 



Using these relations, we can estimate the boosted 
Lorentz factor "/bl- Inside the rarefaction region, the 
gas pressure decreases to that of the contact discontinu- 
ity {pg^n)- Since Pg^D ^ Pg.Ry we immediately obtain the 
upper bound of 7sl, 



7jet( ) 

^Pg,D ' 



(r-i)/r 



< 



^Pg.R> 



1/4 



(13) 



It is interesting to see that 7bl is controlled by the ex- 
ternal pressure Pt,R- The over- pressured jet pushes the 
discontinuity outward, and the external pressure termi- 
nates the boost by stopping the further development of 
the rarefaction structure. The external pressure does no 
mechanical work on the jet fluid. 

Note that the boost does not operate when the jet- 
side pressure becomes nonrelativistic (T < 1). We have 
another restriction from Equation |10[ 



IBL < 7jet( — )• 
^ PL ' 



(14) 



This will replace Equation |13[ when the external pressure 
is too low {pg^B. — > 0). 

We also examine the energy equation. Inside the over- 
pressured {j)g ^ p) and relativistically-moving (47|g( ^ 
1) jet, the fluid energy density is 

E = (7^Wf - Pg) - {-ip)^h. (15) 

Substituting Equation [l5] into Equation |3j we obtain the 
same condition as Equation |9] 

^p{dt + V ■ V)(7/i) + dticip) + V • (jpv) 7/1 

^lP§-^hh)^0. (16) 



Equations 15 and 16 tell us that a specific energy density 
(the energy density per the lab- frame gas density) is con- 
served during the fluid convection. This is because the 
total energy flow {•j'^wtv ~ Aj'^pgv) is much larger than 
the work to expand the jet outward (pgf), and because 
the ideal fluid contains no heat transfer in its proper 
frame. 
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4.2. RMHD theory 

Let us consider the effect of the jet-ahgned magnetic 
field, = (0,0,52). After some algebra in Equations 
[2| we find both the z-momentum and the xz component 
of the stress-energy tensor are unchanged from Iwdrod 
namic ones. Therefore we can utilize Equations [8] and 
We further consider flux conservation. 



IP 



~ const. 



Combining this with Equation we obtain 

^^(r-i)/(2-r) _ ^^^^^ 



(17) 



(18) 



When W2 ~ 1 like the boosted rarefaction region, the 
magnetic pressure approximates ~ |^z- From 

Equations [12] and [18] we construct the pressure condi- 
tion across the tangential discontinuity. 



r2 ^ 2(2-r) 

2 \1bl) 



' Pt.D > Pt,R\ 



r(19) 



The power indexes are both 4 when F = 4/3. Therefore 
we obtain a generalized upper bound, 

1/4 

, ) ■ (20) 



IBL < Ijet 



13 



Note that the total pressure pt replaces the gas pressure 
Pg in Equation 
In the case o: 
(0,B„0) 
tions (e.g 



The perpendicular magnetic field, 
the initial choice of v,, 

bn = K. = 0), »V.^C.^O^ 

main zero (Romero et al. 20051. In this case, the boost 



— simplifies the equa- 
because both v,. and B, re- 



comes from the temporal decrease of the total pressure, 
^'^Wt{D / Dt)vz ~ —dtPt- From Equations [l] and [2] we 
can similarly derive the conservation law. 



D / wt 



)-7P;^(7(/^+^)i'.) =0. (21) 



For simplicity, we consider the magnetically dominated 
limit of b'^/p:^ h (or &^ > Apg). In the jet side {v^ ^ 1) 
we expect jb'^/p ^ const. Combining this with the flux 
conservation 



we expect 



7 



&2 
n2 



const. 



const. 



(22) 



(23) 



The condition across the discontinuity leads to an upper 
bound of 7sL, 



Pt,L 



•Pt,D >Pt.R, 



IBL < Ijet ( 



Pt,L 
Pt.R 



1/2 



(24) 



(25) 



Furthermore, from the polytropic law and Equation 22 
we see that the magnetic pressure decays more rapidly 
than the gas pressure. 



&2 oc p2/r 



^pT- 



(26) 



Consequently, the system behaves similarly as the hydro- 
dynamic case once the gas contribution and the magnetic 
contribution become comparable. Therefore, we usually 
expect intermediate results between Equations [20] and 

m 

Among the two RMHD cases, the boost is more sig- 
nifica nt in the perpendic ular case than in the parallel 
case ( jMizuno et al.[[20"08 ). This is because more electro- 
magnetic energy and momentum are available per a gas 
medium — the jet initially contains larger field energy 
^(-B^ + E"^) and carries additional upward momentum in 
a form of Foynting flux {E x B). We also recall that 
the boost process is related to the pressure decrease, and 
that the magnetic pressure preferably works in the per- 
pendicular directions. 
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Figure 3. Anomalous boost {"yBh/ljet) as a function of the 
total pressure {pt.h/Pt.R)- Three models (HI, Ml, and M2) are 
compare d wi th the theories: Equations |13| and |20| (solid line) and 
Equation |25| (doffed line). 



4.3. Numerical Tests 

In order to verify the scaling theory, we carry out se- 
ries of parameter surveys, by controlling the external 
pressure, pt^R = Pg.R (Table jl]). Figure Js] shows the 
boosted Lorentz factors {^bl) m our RMHETsimulations 
as a function of (pt. L/Pt.R.)- Those values are check ed by 
analytic solutions (Giacomazzo & Rezzolla 2006). For 
example, in the reference cases [PtjJPt.R = 10), the 



theory predicts ^Bh/ljet ^ 1-78 (Equations 13 and 201 



and ^Bb/ljet ^ 3.16 (Equation 25), while we obtain 



iBblljet = 1.67 (HI), 1.59 (Ml), and 2.31 (M2) (see 
also Figure [2];). In general, one can see that the scaling 
laws are in excellent agreement with the boost amplitude 
in the HI and Ml series. The Ml cases are slightly af- 
fected by another limitation (e.g.. Equation 14), due to 
the lower initial temperature {pnjJpL) in the jet. In the 
case of the M2 series. Equation |25| works as a looser up- 
per limit. Since the theory is valid when the magnetic 
pressure dominates in the jet, pt^L ~ it is reason- 

able that we obtain intermediate results in these specific 
cases. 

We perform another parameter survey by reducing the 
jet-side density pL (Table [T]). The results are very simi- 
lar. Since pL ^ Pl, we have even better agreement with 
the theory in the Ml series. 

5. RELEVANCE FOR KINETIC MODELS 
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In this section, we examine the problem from the view- 
point of the kinetic theory. For brevity, we assume that 
the gas moves to the -|-z-direction with a speed of /3 = Vz, 
and we set the particle rest mass to m = 1 . Although the 
RHD theory does not assume a specific dist r ibution func- 
tion, a drifting Maxwellian (Jiittner 1911 Synge 19571 
will be the best starting point: 



f{p)dp oc exp 



7(Po - I^Pz) 
T 



dp, 



(27) 



11/2 



IS 



where p is the particle momentum, po = [1 + (p • p 
the particle energy, and 7, (3 are the fluid bulk properties. 

Shown in Figure H] are momentum-space profiles of 
sample distributio n f unctions. Two samples are gener- 
ated by Equation [27j (1) T = 100 and 7 = 7 and (2) 
T = 70 and 7 = lOsuch that they satisfy Equation [lO| 
We intend to mimic (1) the initial condition in the jet 
and (2) the evolved population in the rarefaction region, 
in model HI. The lab- frame density is set to the same. 

The p2:-pi'ofiles (Figure|4]3) are reasonably different due 
to the thermal spread. In contrast, the Pz-profiles (Figure 
|4^) , which significantly extend to the +pz direction, look 
quite similar. In the left side of the Pz-space {pz ^ 0), 
from Equation [27| and Pz ~ —po, the asymptotic slope 
index s of the distribution F{pz) oc e*^^ yields 



7(1 + 
T 



27 
T ' 



(28) 



We see that the population is quite limited in this side, 
when 7 is large. In the right side, the index s will be 



7(1-/^) 
T 



1 



27r 



const. 



(29) 



Therefore the Pz-profile remains similar in this side, even 
when the "fluid" velocity changes. In addition, since the 
right-side population mainly carries the momentum and 
the energy, the two distributions carry nearly the same 
amount of the momentum and the energy density per the 
lab- frame density, as mentioned by Equations |8] and |16[ 
The relative differences are 0.3% in momentum and 0.6% 
in energy, respectively. 



Important implication of Equation 29 is that the typi- 
cal momentum spread is ~ 27r in the -|-pz-direction. 
Recalling the effective boost condition of T 3> 1, we see 
that a thermal umbrella is much bigger (2T times) than 
the bulk Lorentz factor 7 in the relativistic momentum 
space. 

6. DISCUSSION AND SUMMARY 

As shown in Equation [Tj the anomalous bulk boost 
comes from the temporal decrease of relativistic pres- 
sure. From the energy viewpoint, the term transports 
the internal ener gy to that of the bulk m otion {pg 7), 
as mentioned by Aloy & Mimica (20081. The internal- 



to-bulk energy transport is somewhat counter-intuitive, 
however, it is a logical consequence of the relativistic fluid 
formalism. 

The site of the boost is the rarefaction region. The 
rarefaction wave involves the temporal pressure decrease 
behind its wave front and there is a room for the con- 
vective fluid motion (Equation [sl . In contrast, neither 
conditions are satisfied around the shocks. The anoma- 
lous boost does not occur on the other side of the con- 
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Figure 4. (Color online) (a) Ideal gas distribution functions for 
(1) T = 100 and 7 = 7 {solid line) and (2) T = 70 and 7 = 10 
{dotted line) in the p^^-spacc. (b) The same, but in the pz-space. 



tact/tangential discontinuity nor will it occur when an- 
other shock replaces the rarefaction wave. Therefore, 
the transitio n from the shock regime t o the rarefaction 
wave regime ( RezzoUa & Zanottil|2002 1 would be a criti- 



cal condition for the problem. Similar boost in the nor- 
mal direction has recently been reported i n magnetically- 



domi nated rarefaction region as well (Mizuno et al. 



2009| ) 

Xnother explan ation is a relativistic fr ee expansion 
in the jet frame (Komissarov et al. 20091. When the 



relativistically strong pressure pushes the gas outward 
against the external medium, the lateral expansion can 
be relativistic in the jet frame. Then, the Lorentz fac- 
tor in the observer frame yields ^bl ~ 7iet(l — f'^)"^/^, 
where v' is the expansion speed in the jet frame. We 
expect that the term —v'dt'Pt enhances such expansion 
in the rarefaction region, and that the relevant boost is 
projected into the tangential boost in the observer frame. 
Strictly speaking, a ID problem in the observer frame is 
no longer identical to that in the jet frame, because a ID 
expansion of the discontinuity front in the -|-a;-direction 
is projected to the oblique direction in the jet frame. The 
two problems start differently and therefore the situation 
is more complicated. 

A potential limitation is that multi-dimensional insta- 
bilities may modulate the ID evolution. Especially, the 
relativistic Kelvin-Helmholtz (KH) instabilities will be 
relevant. In th e regime of our interest, t he increasing 
Lorentz factor fTurland & Scheuer 1976 Blandford &j 
Pringle_1976; Bodo et jil. 2004) and the how-aligned mag- 
netic held fOsmanov et a l. 2008) suppres s the K H mode; 
for instance, it we employ Bodo et al. (2004rs stabil- 
ity condition of jjet > (1 -I- 2 cos~^ 9) in our RHD jet 
iljet = 7), where 6 is the angle between the jet flow 
and the wavevector, the instability is allowed only in 
the quasi-transverse direction. On the other hand, shear 
layers with density asymmetry are known to be substan- 
tially KH-unstable. Once the KH vortex develops, the 
subsequent turbulence is likely to smooth the sharp lat- 
eral structure. While ID-like sign atures have been f ound 
in some three-dimensional RHD (Aloy et al."2005) and 
two-dimensional RM HD simu lations (Mizuno et al . 20081 
Tchekhovskoy et al. 2009; K omissarov et al.|2U09 ), "inter- 



ference with the KH and other instabilities needs further 
investigation. 

In addition, we need to keep in mind that the entire 
process depends on the ideal fluid assumption. In or- 
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der to justify it, coUisional or other scattering processes 
have to relax the gas much quicker than the dynamical 
timescale. However, those are difficult conditions espe- 
cially in the jet side, where the physical processes look 
even slower by the rclativistic effect. In Section 5, we 
show that the fluid bulk speed is considerably smaller 
than a wide thermal spread in the momentum profile, 
when the boost operates. We think that the counter- 
intuitive force may be just enforced by the ideal fiuid 
assumption: i.e. the anomalous fluid acceleration may 
be an artifact of an expedient isotropic fluid velocity. In 
the real world, we expect that non-ideal effects such as 
the heat flow play roles. In fact, the system involves large 
gradient of the pressure and the temperature in the rar- 
efaction regions and around the discontinuities. In the 
high-temperature regime of T 3> 1, the energy and mo- 
mentum balances are mainly controlled by the pressure 
parts (the internal energy or the enthalpy flux), which 
can be sensitive to the local gas distribution functions. 

In sum mary, we examined the ID anomalous relativ is- 
tic boost ( |Aloy fc Rezzona||2006[ [Mizuno et aT1|2008[ ) at 
the lateral boundary of relativistic jets. We numerically 
and theoretically conflrmed that the anomalous boost oc- 
curs in the RHD and RMHD regimes. We further derived 
simple scaling laws for the accelerated Lorentz factor. 



IBL < Ijet 



1/4 (hydro, parallel) 
s = 1/2 (perpendicular) 



We also note that the process operates in an ideal fluid. 
The non-ideal effects (heat flow etc.) as well as multi- 
dimensional effects are left for future works. We hope 
that this work will be a basic piece for the boundary 
problems in rclativistic jets and the relevant simulations. 
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